Evidence for hybridization-driven heteroplasmy maintained across generations in a ricefish endemic to a Wallacean ancient lake.

Heteroplasmy, the presence of multiple mitochondrial DNA (mtDNA) haplotypes within cells of an individual, is caused by mutation or paternal leakage. However, heteroplasmy is usually resolved to homoplasmy within a few generations because of germ-line bottlenecks; therefore, instances of heteroplasmy are limited in nature. Here, we report heteroplasmy in the ricefish species Oryzias matanensis, endemic to Lake Matano, an ancient lake in Sulawesi Island, in which one individual was known to have many heterozygous sites in the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene. In this study, we cloned the ND2 gene for some additional individuals with heterozygous sites and demonstrated that they are truly heteroplasmic. Phylogenetic analysis revealed that the extra haplotype within the heteroplasmic O. matanensis individuals clustered with haplotypes of O. marmoratus, a congeneric species inhabiting adjacent lakes. This indicated that the heteroplasmy originated from paternal leakage due to interspecific hybridization. The extra haplotype was unique and contained two non-synonymous substitutions. These findings demonstrate that this hybridization-driven heteroplasmy was maintained across generations for a long time to the extent that the extra mitochondria evolved within the new host.

Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples.

BACKGROUND: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. RESULTS: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. CONCLUSIONS: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

Quantifying mitochondrial heteroplasmy diversity: A computational approach.

Biodiversity plays a pivotal role in sustaining ecosystem processes, encompassing diverse biological species, genetic types and the intricacies of ecosystem composition. However, the precise definition of biodiversity at the individual level remains a challenging endeavour. Hill numbers, derived from Rényi's entropy, have emerged as a popular measure of diversity, with a recent unified framework extending their application across various levels, from genetics to ecosystems. In this study, we employ a computational approach to exploring the diversity of mitochondrial heteroplasmy using real-world data. By adopting Hill numbers with q = 2, we demonstrate the feasibility of quantifying mitochondrial heteroplasmy diversity within and between individuals and populations. Furthermore, we investigate the alpha diversity of mitochondrial heteroplasmy among different species, revealing heterogeneity at multiple levels, including mitogenome components and protein-coding genes (PCGs). Our analysis explores large-scale mitochondrial heteroplasmy data in humans, examining the relationship between alpha diversity at the mitogenome components and PCGs level. Notably, we do not find a significant correlation between these two levels. Additionally, we observe significant correlations in alpha diversity between mothers and children in blood samples, exceeding the reported R2 value for allele frequency correlations. Moreover, our investigation of beta diversity and local overlay similarity demonstrates that heteroplasmy variant distributions in different tissues of children more closely resemble those of their mothers. Through systematic quantification and analysis of mitochondrial heteroplasmy diversity, this study enhances our understanding of heterogeneity at multiple levels, from individuals to populations, providing new insights into this fundamental dimension of biodiversity.

Whole mitogenome sequencing uncovers a relation between mitochondrial heteroplasmy and leprosy severity.

BACKGROUND: In recent years, the mitochondria/immune system interaction has been proposed, so that variants of mitochondrial genome and levels of heteroplasmy might deregulate important metabolic processes in fighting infections, such as leprosy. METHODS: We sequenced the whole mitochondrial genome to investigate variants and heteroplasmy levels, considering patients with different clinical forms of leprosy and household contacts. After sequencing, a specific pipeline was used for preparation and bioinformatics analysis to select heteroplasmic variants. RESULTS: We found 116 variants in at least two of the subtypes of the case group (Borderline Tuberculoid, Borderline Lepromatous, Lepromatous), suggesting a possible clinical significance to these variants. Notably, 15 variants were exclusively found in these three clinical forms, of which five variants stand out for being missense (m.3791T > C in MT-ND1, m.5317C > A in MT-ND2, m.8545G > A in MT-ATP8, m.9044T > C in MT-ATP6 and m.15837T > C in MT-CYB). In addition, we found 26 variants shared only by leprosy poles, of which two are characterized as missense (m.4248T > C in MT-ND1 and m.8027G > A in MT-CO2). CONCLUSION: We found a significant number of variants and heteroplasmy levels in the leprosy patients from our cohort, as well as six genes that may influence leprosy susceptibility, suggesting for the first time that the mitogenome might be involved with the leprosy process, distinction of clinical forms and severity. Thus, future studies are needed to help understand the genetic consequences of these variants.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Vanadium Pentoxide Exposure Causes Strain-Dependent Changes in Mitochondrial DNA Heteroplasmy, Copy Number, and Lesions, but Not Nuclear DNA Lesions.

Interstitial lung diseases (ILDs) are lethal lung diseases characterized by pulmonary inflammation and progressive lung interstitial scarring. We previously developed a mouse model of ILD using vanadium pentoxide (V2O5) and identified several gene candidates on chromosome 4 associated with pulmonary fibrosis. While these data indicated a significant genetic contribution to ILD susceptibility, they did not include any potential associations and interactions with the mitochondrial genome that might influence disease risk. To conduct this pilot work, we selected the two divergent strains we previously categorized as V2O5-resistant C57BL6J (B6) and -responsive DBA/2J (D2) and compared their mitochondrial genome characteristics, including DNA variants, heteroplasmy, lesions, and copy numbers at 14- and 112-days post-exposure. While we did not find changes in the mitochondrial genome at 14 days post-exposure, at 112 days, we found that the responsive D2 strain exhibited significantly fewer mtDNA copies and more lesions than control animals. Alongside these findings, mtDNA heteroplasmy frequency decreased. These data suggest that mice previously shown to exhibit increased susceptibility to pulmonary fibrosis and inflammation sustain damage to the mitochondrial genome that is evident at 112 days post-V2O5 exposure.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Negative correlation between organ heteroplasmy, particularly hepatic heteroplasmy, and age at death revealed by post-mortem studies of m.3243A > G cases.

Mitochondrial DNA m.3243A > G mutation causes mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and its associated multi-organ disorders, including diabetes. To clarify associations between m.3243A > G organ heteroplasmy and clinical phenotypes, including the age at death, we combined genetic and pathological examinations from seven unreported and 36 literature cases of autopsied subjects. Clinical characteristics of subjects were as follows: male, 13; female, 28; unknown, 2; the age at death, 36.9 ± 20.2 [4-82] years; BMI, 16.0 ± 2.9 [13.0-22.3]; diabetes, N = 21 (49%), diabetes onset age 38.6 ± 14.2 years; deafness, N = 27 (63%); stroke-like episodes (StLEp), N = 25 (58%); congestive heart failure (CHF), N = 15 (35%); CHF onset age, 51.3 ± 14.5 years. Causes of death (N = 32) were as follows: cardiac, N = 13 (41%); infection, N = 8 (25%); StLEp, N = 4 (13%); gastrointestinal, N = 4 (13%); renal, N = 2 (6%); hepatic, N = 1 (2%). High and low heteroplasmies were confirmed in non-regenerative and regenerative organs, respectively. Heteroplasmy of the liver, spleen, leukocytes, and kidney for all subjects was significantly associated with the age at death. Furthermore, the age at death was related to juvenile-onset (any m.3243A > G-related symptoms appeared before 20) and stroke-like episodes. Multiple linear regression analysis with the age at death as an objective variable showed the significant contribution of liver heteroplasty and juvenile-onset to the age at death. m.3243A > G organ heteroplasmy levels, particularly hepatic heteroplasmy, are significantly associated with the age at death in deceased cases.

PHFinder: assisted detection of point heteroplasmy in Sanger sequencing chromatograms.

Heteroplasmy is the presence of two or more organellar genomes (mitochondrial or plastid DNA) in an organism, tissue, cell or organelle. Heteroplasmy can be detected by visual inspection of Sanger sequencing chromatograms, where it appears as multiple peaks of fluorescence at a single nucleotide position. Visual inspection of chromatograms is both consuming and highly subjective, as heteroplasmy is difficult to differentiate from background noise. Few software solutions are available to automate the detection of point heteroplasmies, and those that are available are typically proprietary, lack customization or are unsuitable for automated heteroplasmy assessment in large datasets. Here, we present PHFinder, a Python-based, open-source tool to assist in the detection of point heteroplasmies in large numbers of Sanger chromatograms. PHFinder automatically identifies point heteroplasmies directly from the chromatogram trace data. The program was tested with Sanger sequencing data from 100 humpback whales (Megaptera novaeangliae) tissue samples with known heteroplasmies. PHFinder detected most (90%) of the known heteroplasmies thereby greatly reducing the amount of visual inspection required. PHFinder is flexible and enables explicit specification of key parameters to infer double peaks (i.e., heteroplasmies).

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans.

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

CmVCall: An automated and adjustable nanopore analysis pipeline for heteroplasmy detection of the control region in human mitochondrial genome.

Genetic associations between human mitochondrial DNA (mtDNA) heteroplasmy and mitochondrial diseases, aging, and cancer have been elaborated, contributing a lot to the further understanding of mtDNA polymorphic spectrum in anthropology, population, and forensic genetics. In the past decade, heteroplasmy detection using Sanger sequencing and next generation sequencing (NGS) was hampered by the former's inefficiency and the latter's inherent bias due to amplification and mapping of short reads, respectively. Nanopore sequencing stands out for its ability to yield long contiguous segments of DNA, providing a new insight into heterogeneity authentication. In addition to MinION from Oxford Nanopore Technologies, an alternative nanopore sequencer QNome (Qitan Technology) has also been applied to various biological research and the forensic applicability of this platform has been proved recently. In this study, we evaluated the performance of four commonly used variant callers in the heterogeneity authentication of the control region of human mtDNA based on simulations of different ratios generated by mixing QNome nanopore sequencing reads of two synthetic sequences. Then, an open-source and python-based nanopore analytics pipeline, CmVCall was developed and incorporated multiple programs including reads filtering, removal of nuclear mitochondrial sequences (NUMTs), alignment, optional 'Correction' mode, and heterogeneity identification. CmVCall can achieve high precision, accuracy, and recall of 100%, 99.9%, and 92.3% with a 5% heteroplasmy level in 'Correction' mode. Moreover, blood, saliva, and hair shaft samples from monozygotic (MZ) twins were used for heterogeneity evaluation and comparison with the NGS data. Results of MZ twin samples showed that CmVCall could identify more point heteroplasmy sites, revealing significant levels of inter- and intra-individual mtDNA polymorphism. In conclusion, we believe that this analysis pipeline will lay a solid foundation for the development of a comprehensive nanopore analysis pipeline targeting the whole mitochondrial genome.

Genetic diversity of thrips populations on Allium species around the world.

Thrips are a serious pest in many crops. In onion cultivation, Thrips tabaci is the most important, but not the only thrips species causing damage. We investigated which thrips species affects onion and related species worldwide, how much genetic variation there is within T. tabaci populations, and how this evolves. Furthermore, we determined the reproductive mode and the correlation between the genetic and geographic distances. Thrips samples from infested onions or related species were obtained from 14 different locations worldwide. Species and haplotypes were determined through DNA barcoding with the mitochondrial Cytochrome Oxidase subunit I (COI) gene. Thrips tabaci was the most commonly observed species, but Scirtothrips dorsalis, Thrips palmi, Frankliniella intonsa, Frankliniella occidentalis and Frankliniella tenuicornis were also found, especially at the beginning of the growing seasons and depending on the location. The Nei's genetic distance within T. tabaci was less than 5% and the haplotypes were clustered into two phylogenetic groups, each linked to a specific mode of reproduction, thelytokous or arrhenotokous. Thelytokous thrips were more common and more widely distributed than arrhenotokous thrips. A high percentage of heteroplasmy was detected in the arrhenotokous group. Heteroplasmic thrips were only found in populations where thelytokous and arrhenotokous were present in sympatry. Some T. tabaci haplotypes were present in high frequency at several sampled locations. No correlation was found between the genetic and geographic distances, which points to anthropic activities spreading thrips haplotypes throughout the world.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.).

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

The complete mitochondrial genome of Trigonisca nataliae (Hymenoptera, Apidae) assemblage reveals heteroplasmy in the control region.

The evolution of mitochondrial genomes in the stingless bees is surprisingly dynamic, making them a model system to understand mitogenome structure, function, and evolution. Out of the seven mitogenomes available in this group, five exhibit atypical characteristics, including extreme rearrangements, rapid evolution and complete mitogenome duplication. To further explore the mitogenome diversity in these bees, we utilized isolated mtDNA and Illumina sequencing to assemble the complete mitogenome of Trigonisca nataliae, a species found in Northern Brazil. The mitogenome of T. nataliae was highly conserved in gene content and structure when compared to Melipona species but diverged in the control region (CR). Using PCR amplification, cloning and Sanger sequencing, six different CR haplotypes, varying in size and content, were recovery. These findings indicate that heteroplasmy, where different mitochondrial haplotypes coexist within individuals, occurs in T. nataliae. Consequently, we argue that heteroplasmy might indeed be a common phenomenon in bees that could be associated with variations in mitogenome size and challenges encountered during the assembly process.

Modeling of mitochondrial genetic polymorphisms reveals induction of heteroplasmy by pleiotropic disease locus 10398A>G.

Mitochondrial (MT) dysfunction has been associated with several neurodegenerative diseases including Alzheimer's disease (AD). While MT-copy number differences have been implicated in AD, the effect of MT heteroplasmy on AD has not been well characterized. Here, we analyzed over 1800 whole genome sequencing data from four AD cohorts in seven different tissue types to determine the extent of MT heteroplasmy present. While MT heteroplasmy was present throughout the entire MT genome for blood samples, we detected MT heteroplasmy only within the MT control region for brain samples. We observed that an MT variant 10398A>G (rs2853826) was significantly associated with overall MT heteroplasmy in brain tissue while also being linked with the largest number of distinct disease phenotypes of all annotated MT variants in MitoMap. Using gene-expression data from our brain samples, our modeling discovered several gene networks involved in mitochondrial respiratory chain and Complex I function associated with 10398A>G. The variant was also found to be an expression quantitative trait loci (eQTL) for the gene MT-ND3. We further characterized the effect of 10398A>G by phenotyping a population of lymphoblastoid cell-lines (LCLs) with and without the variant allele. Examination of RNA sequence data from these LCLs reveal that 10398A>G was an eQTL for MT-ND4. We also observed in LCLs that 10398A>G was significantly associated with overall MT heteroplasmy within the MT control region, confirming the initial findings observed in post-mortem brain tissue. These results provide novel evidence linking MT SNPs with MT heteroplasmy and open novel avenues for the investigation of pathomechanisms that are driven by this pleiotropic disease associated loci.

Effects of Parkin on the Mitochondrial Genome in the Heart and Brain of Mitochondrial Mutator Mice.

Mitochondrial dysfunction has been implicated in neurodegenerative diseases like Parkinson's disease (PD). This study investigates the role of Parkin, a protein involved in mitochondrial quality control, and strongly linked to PD, in the context of mitochondrial DNA (mtDNA) mutations. Mitochondrial mutator mice (PolgD257A/D257A ) (Polg) are used and bred with Parkin knockout (PKO) mice or mice with disinhibited Parkin (W402A). In the brain, mtDNA mutations are analyzed in synaptosomes, presynaptic neuronal terminals, which are far from neuronal soma, which likely renders mitochondria there more vulnerable compared with brain homogenate. Surprisingly, PKO results in reduced mtDNA mutations in the brain but increased control region multimer (CRM) in synaptosomes. In the heart, both PKO and W402A lead to increased mutations, with W402A showing more mutations in the heart than PKO. Computational analysis reveals many of these mutations are deleterious. These findings suggest that Parkin plays a tissue-dependent role in regulating mtDNA damage response, with differential effects in the brain and heart. Understanding the specific role of Parkin in different tissues may provide insights into the underlying mechanisms of PD and potential therapeutic strategies. Further investigation into these pathways can enhance the understanding of neurodegenerative diseases associated with mitochondrial dysfunction.

Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease.

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multisystem disease. The wide range of manifestations observed in patients with mitochondrial disease results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify nonrandom distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA-Seq, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profiled the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and from healthy control donors. Utilizing the retina as a model for complex multilineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm-derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm-derived lineage, namely the vasculature of the choroid, was near homoplasmic for the WT allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the nonrandom nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

Discrimination of monozygotic twins using mtDNA heteroplasmy through probe capture enrichment and massively parallel sequencing.

Differentiating between monozygotic (MZ) twins remains difficult because they have the same genetic makeup. Applying the traditional STR genotyping approach cannot differentiate one from the other. Heteroplasmy refers to the presence of two or more different mtDNA copies within a single cell and this phenomenon is common in humans. The levels of heteroplasmy cannot change dramatically during transmission in the female germ line but increase or decrease during germ-line transmission and in somatic tissues during life. As massively parallel sequencing (MPS) technology has advanced, it has shown the extraordinary quantity of mtDNA heteroplasmy in humans. In this study, a probe hybridization technique was used to obtain mtDNA and then MPS was performed with an average sequencing depth of above 4000. The results showed us that all ten pairs of MZ twins were clearly differentiated with the minor heteroplasmy threshold at 1.0%, 0.5%, and 0.1%, respectively. Finally, we used a probe that targeted mtDNA to boost sequencing depth without interfering with nuclear DNA and this technique can be used in forensic genetics to differentiate the MZ twins.

Mitochondrial Heteroplasmy and PCR Amplification Bias Lead to Wrong Species Delimitation with High Confidence in the South American and Antarctic Marine Bivalve Aequiyoldia eightsii Species Complex.

The species delimitation of the marine bivalve species complex Aequiyoldia eightsii in South America and Antarctica is complicated by mitochondrial heteroplasmy and amplification bias in molecular barcoding. In this study, we compare different data sources (mitochondrial cytochrome c oxidase subunit I (COI) sequences; nuclear and mitochondrial SNPs). Whilst all the data suggest that populations on either side of the Drake Passage belong to different species, the picture is less clear within Antarctic populations, which harbor three distinct mitochondrial lineages (p-dist ≈ 6%) that coexist in populations and in a subset of individuals with heteroplasmy. Standard barcoding procedures lead to amplification bias favoring either haplotype unpredictably and thus overestimate the species richness with high confidence. However, nuclear SNPs show no differentiation akin to the trans-Drake comparison, suggesting that the Antarctic populations represent a single species. Their distinct haplotypes likely evolved during periods of temporary allopatry, whereas recombination eroded similar differentiation patterns in the nuclear genome after secondary contact. Our study highlights the importance of using multiple data sources and careful quality control measures to avoid bias and increase the accuracy of molecular species delimitation. We recommend an active search for mitochondrial heteroplasmy and haplotype-specific primers for amplification in DNA-barcoding studies.